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Abstract 

We compute the cross section for XX ~^ QQ9 a t order o? s jM~ arising from interference between the 
tree-level and loop-induced processes. This interference term is the same order in a s as XX ~> 99\ 
for mass degenerate squarks M~ = M~ = M~ we find v xe \a- m i = [— 2m?,/3M~] v Ie \a(xx ~^ 99)- 
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I. INTRODUCTION 



The presence of non-baryonic dark matter in the universe is compelling evidence for 
physics beyond the Standard Model. Supersymmetry (SUSY) provides an especially attrac- 
tive explanation with the lightest supersymmetric particle (LSP), usually a neutralino x, as 
the dark matter candidate; for a recent review, see Ref. Q]. 

We focus in this paper on the QCD corrections to neutralino annihilation. We therefore 
consider only the annihilation processes involving quarks and/or gluons in the final state. 
We also assume that the LSP is largely gaugino, as is often the case in mSUGRA models 
in this case neutralino annihilation via s-channel Higgs or Z exchange is suppressed, since 
these particles require a Higgsino admixture to couple to the LSP. We thus focus on processes 
that involve an internal squark exchange, as shown in Fig. ^ 

The leading contribution to neutralino annihilation via exchange of a squark of mass 
M~, shown in Fig. ^a), can be reduced to an effective vertex described by a dimension-six 
operator suppressed by M~, 

C = (c/Mj)0 6 , O e = (n,75X)(?7"75?), (1) 

where the dimensionless constant c contains the relevant couplings of the process. If the 
annihilation occurs at rest, the neutralino spinors reduce to a spin-singlet combination, 
X7a*75X — > K^xilsl v2> m x)Xi where K is the center-of-mass momentum of the two-neutralino 
system. This result can be understood by considering that the Majorana nature of neutrali- 
nos makes the initial state behave as an effective pseudoscalar (details are given in Sec. ITTj) . 
The operator Oq can then be written as the divergence of the axial- vector current, 

[d,(qYl 5 q)] ■ (2) 

In the massless quark limit, m q = 0, the axial vector current is conserved at tree level, 
^(97 M 75?) = 0, an d all tree-level dimension-six amplitudes vanish in this limit. This is the 
well known partially-conserved axial current (PCAC) condition. 
The quark mass suppression can be lifted in two ways: 

1. by going beyond leading order in a s to include the correction to the dimension-six 
operator that involves the anomalous triangle diagram; 

2. by going to dimension-eight by including hard gluon radiation. 
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FIG. 1: Feynman diagrams that contribute to neutralino annihilation to quarks /gluons. Shown are 
the tree- level diagram for XX ~^ QQ ( a )> typical one- loop diagrams for XX ~~ > 99 (b) and XX ~^ QQ9 
(c), and the three tree-level diagrams contributing to the dimension-eight operator for XX ~> QQ9 
(d-f). 



The anomalous triangle diagram is shown in Fig. E^b) (the quark triangle appears at 
dimension-six when the heavy squark line is shrunk to a point). This diagram contributes 
t° XX ~^ 99 an d yields a cross section parametrically of order a 2 s m 2 /M~. The analogous 



process XX ~^ 77 was nrs t studied using the anomaly equation in Refs. I3j, m; a modern 
computation of XX 99 i n the v re \ = limit with arbitrary neutralino composition and 
squark mixing, and keeping m q / ^ in the loop, was performed in Ref. Q. Ref. Q extended 
the calculation to arbitrary t> rol , with results presented in a set of numerical codes 

The next-to-leading order (NLO) QCD corrections to XX 99 were recently computed 
in Ref. jsl in the approximation m q > = and M~ ^> m x . These corrections are of order 
alm 2 /M~, and provide about a 60% enhancement over the leading order (LO) XX 99 



to a dimension-eight 



cross section for typical parameter values 

The tree-level process XX ~^ QQ9 corresponds in the limit m q 



operator, with the amplitude suppressed by 1/M~. The relevant diagrams are shown in 
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Process Order Ref. Approximation 

XX ~^ 99-, Dini-6 a 2 s m^/M~ [5] exact for v re \ = 

XX ~~ * QQ9> Dim-8 a s m x /M~ [5] exact for m q = 

XX ~^ QQ9i Dim-6/8 interference a 2 s m x /M~ [this paper] leading 1/M~ term, m q ^ q i 



XX ~~ * 99 a t NLO, Dim-6 o? s m \/M~ [8] leading 1/M~ term, m q / = 



TABLE I: Relevant neutralino annihilation calculations and their relative suppression by powers 
of a s and vntjM-. 



Fig. md-f). The similar process XX ffl was calculated in Refs. j4j, |9( for a photino 
LSP. An exact tree-level calculation of XX ~~ * QQ9 with arbitrary neutralino composition and 
squark mixing, for m q = and v re \ = 0, was done in Ref. and has a resulting cross section 
of order a s m x /M~. 

In this paper we compute the interference term between the tree-level XX QQ9 diagrams 
of Figs. HKd,e,f) an d the one-loop XX QQ9 process of Fig. HJ^c). This interference term 
is parametrically of order a 2 s m x /M~. Counting a s and m x /M~ as comparable suppression 
factors, this interference term is "the same order" as both the dimension-eight process and 
the NLO QCD corrections to XX ~^ 99) as summarized in Table H] Counting only the order 
in a s , the interference term is the same order as XX 99- 

Throughout this paper we assume that the LSP is lighter than the top quark, so that 



XX ~ > is kinematically inaccessible. (T 
threshold were recently computed in Ref. 



le QCD corrections to XX Z*,h* — > tt near 



id].) We take the five light quark species to be 



massless. As noted earlier, we also assume that the neutralinos annihilate at rest and that 
the LSP is gaugino-like (the coupling of the Higgsino component to quark-squark pairs is 
proportional to the quark mass, taken here to be zero). We also assume that the squarks 
are heavy compared to the LSP, M~ ^> m x , and work to leading order in the expansion in 
m x /M~. In the calculation of the quark/squark loop in Fig. ^c) we sum over five mass- 
less internal quark species, and ignore the top quark contribution; the massive quark loop 
decouples very quickly with decreasing m x /m t . 
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II. CALCULATION 



In the limit of zero relative neutralino velocity, v re \ = 0, a neutralino pair behaves as a 
pseudoscalar due to the Majorana nature of the particles. In particular, the antisymmetrized 
wave function of the initial-state identical particles xXi °f total momentum K, can be 
reduced to a projection operator for each of the four spin combinations of the initial-state 
neutralinos: 

Ul (K/2)v 2 (K/2) - u 2 {K/2)v 1 {K/2) = (m x + #/2) 75 for | ->| <-) or - | <-)| ->, 
u 1 (K/2)v 2 (K/2)-u 2 (K/2)v 1 (K/2) = for | ->>| or | <-)]<-). (3) 

The cross section is then obtained by averaging over the squares of the amplitudes of the 
two contributing spin states. 

One can alternatively use the spin singlet initial state, 

(I ->>!«-> -I «->!-» />/2. (4) 

which produces an additional v2 in the antisymmetrized spinor combination and in the 
resulting matrix elements: 

Ul (K/2)v 2 (K/2) - u 2 (K/2)v 1 (K/2) = V2(m x + ft/ 2) 7s (spin singlet). (5) 

Either approach leads to the same overall result for the annihilation cross section; we use 
Eq. JSJ) with the additional \/2 in what follows in order to obtain the matrix element for the 
spin-singlet state. (After squaring the matrix element, one must still average over the four 
initial neutralino spin states, yielding a factor of 1/4 in the cross section.) 

A. Dimension-six amplitude 

In the m q = limit, the leading contribution to the dimension-six amplitude arises at 
order ot 2 s from the loop-induced process XX ~~ > 99- I n the language of the partially conserved 
axial current, this can be expressed as the anomaly equation 

d,(qYl 5 q) = 2m q qi l5 q + ^-G$& a ^, (6) 

where \G ilv = e^ ua pG al3 is the dual of the color field strength tensor. For m q «C m x , the 
first term on the right-hand side of Eq. (jUJ) can be neglected; the dimension-six amplitude 
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can then be expressed in the form 



M loop (xx - 99) = (f^H &™)^G$Gtor 



(7) 



The Feynman diagrams for \X ~ * 99 consist of Fig. IHb), the analogous diagram with one 
neutralino and one gluon line interchanged, and the corresponding diagrams with the gluons 
or the neutralinos crossed. The amplitude for \X ~ * 99 i n the massless quark limit is given 
by 



M 



loop 



0£ s x - 



\9r 



1r 



+ 



\9e\ 
Mi, 

1l 



16 



livxxfi *a 



(ki)el a (k 2 )k la k 2 p, 



(8) 



where the sum runs over the five light (massless) quark flavors q'. An expression for the 
expansion of A^i oop in powers of rri l/M~ is given in the appendix. The xq'q'rl couplings g r / 



for right and left quark helicities, respectively, are 

g r = -V2N u g'Q, g £ = -V2N n gf(T 3 - Q) + V2N 12 gT 3 . 



(9) 



Here T3 is the quark isospin, Q is the quark electric charge, g and g' are the weak couplings 
under SU (2) L and Z7(l)y, respectively, and N u and Ni 2 are the bino and wino components of 
the neutralino as defined in Ref. 3]. The matrix element in Eq. (jSJ) results in an annihilation 
cross section for XX ~^ 99 °f 



^relO"dim6 



327T 3 X 



E 



\9r 



Ml, 



+ 



\9tr 

1 L 



(10) 



For massless quarks in the loop, taking one of the final-state gluons off-shell by an amount 
<? 2 7^ yields a relative shift in the loop integral by q 2 /M~; in particular, this shift is higher 
order in the 1/M~ expansion than the leading term. To leading 1/M~ order it is then 
straightforward to extend the loop amplitude to include one off-shell gluon splitting into a 
quark- ant iquark pair, shown schematically in Fig. dc). The matrix element is 



M 



split 



9s 
V2 



E 



\9r 



Ml, 

1r 



+ 



\9tT 
Ml, 



le 



fiuafi *c 



2n( qi + q 2 y 

where qi, q 2 , and ^3 are the momenta of the final-state quark, antiquark, and gluon, respec- 
tively, and the sum runs over the five light (massless) quarks in the loop. Note that the 
gluon splitting couples the box diagram to a vectorlike quark current. 

The genuine anomaly loop is a closed fermion loop which has a negative sign. The loop in 
the present case is basically a box with one side as the scalar quark line, so it is not a closed 



fermion loop. However, it does have a relative negative sign with respect to the dimension- 
eight amplitude (see the next subsection) because the two fermion lines are topologically 
twisted - replacing the neutralino spinor pair with the zero-velocity result from Eq. (J3J) 
requires contracting the two spinors together, yielding an extra overall minus sign. The 
situation is similar to that of the two amplitudes in Bhabha scattering. 

Using Eq. we can compute the zero-velocity neutralino annihilation cross section 
from the gluon-splitting process. Squaring the matrix element, summing over final-state 
polarizations and colors and averaging over initial-state polarizations, we obtain 

2 



pols 



|2 |„2| 

I + \9t\ 



Mi, Mi, 
. 1r 1l 



qi ■ <?2 



for a single (massless) final-state quark flavor. 

Note the divergence in Eq. (fT2J) as the gluon propagator goes on shell, (gi + q 2 ) 2 = 
• q 2 ~ 0. This is the usual divergence that appears in NLO calculations when a final- 
state parton is soft or collinear; in this case the divergence is cut off by the quark mass and 

n 

gives rise to a \og{m 2 c / m 2 ^ enhanced term in the total cross section |9[. This logarithmic 
term is precisely canceled by the renormalization of the strong coupling due to the quark 
bubble that appears in the virtual part of the NLO correction to XX ~^ 99 jS| • This is the 
familiar cancellation of logarithmic divergences guaranteed by the Kinoshita-Lee-Nauenberg 
theorem jl^ . 



B. Dimension-eight amplitude 

The dimension-eight amplitude comes from computing the tree-level process XX ~^ QQ9 
shown in Fig. ^d-f ) in the massless quark limit, m q = 0. Because the XX ~^ <?<? amplitude 
[Fig. W[&)] is zero in this limit, the gluon radiation diagram contains no soft or collinear 
divergences. The tree-level matrix element has the general form, 

iM tTee = (i) 5 g s \g r AMqi)T c A c P R , L v(q 2 ), (13) 

where the amplitude A c is given below. The i in front of .Mtree follows the usual Feynman 
rule convention. The factor (i) 5 counts three vertices and two propagators. The SU(3) 
matrix T c links the outgoing gluon with color label c (contained inside A c ) to the qq system. 
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The chirality projection operators Pr : l = (1 ±7s)/2 reflect the fact that the squarks couple 
to quarks of specific helicity. The pieces of the amplitude are A c = A\ + A 2 + with 

1 K 



A\ 



4\ + 93 V2 



75 



(K/2 - q 2 y - Mf R L 



<1R,L 



(K/2 - q 2 y - M| 



-K/2 + qi y-M- 



QR,L 



-K/2 + qi y-M- 



QR,L 



(qi - Q2) 



-k/2 + qi y 



Ml 

1R,L 



(K/2 - q 2 y 



Ml 

QR,L 



(14) 



where A\, A 2 , and A\ correspond to cases that the gluon is radiated from the quark line, 
the antiquark line and the squark line, as in Figs. ^d,e,f), respectively. The outgoing four- 
momenta of the quark, antiquark and gluon are denoted qi, q2 and #3, respectively, and K is 
the incoming four-momentum of the initial two-neutralino system. The Dirac equation was 
applied to A\ and A 2 in the massless quark limit. Combining these three contributions, we 
get 



A c 



1 



(q 2 - qi ).K<f c + ( qi -q 2 )-e* c K 



V2 



[-K/2 + qi y-M^ L 



(K/2 - q 2 f - M- 



T75- 



Qr,l 



(15) 



Note that the pieces with only one squark propagator in the denominator have canceled in 
the massless quark limit. The resulting matrix element is 

(q2~qi)-q3(f c + (qi-q2)-e* c (/ 3 



M 



tree 



9s I 1 2 _/ \rrc 

'~^\9r,e\ u(qi)T 



-K/2 + qi y 



M2 

Qr,l 



(K/2 - q 2 y 



Mi 

<1R,L 



l5PR,Lv(q 2 ), 



(16) 



for the diagrams involving q^ and qi, respectively. 
Using the Chisholm identity, 



0°V - 9^ 



7 a 7V-/ M 7 a -^ Q/3 7,75, 



with 



,0123 



+ 1, 



(17) 



together with the Dirac equation and the polarization condition e* = for the external 
gluon, we obtain an amplitude which has a similar tensor structure to that of the dimension- 
six amplitude in Eq. (|TT|) : 

ie^e*;q 3a ( qi + q 2 ) p 



Mtrcc = \g r ,e\ u(qi)T c -f u P R)L v(q 



-K/2 + gi y 



Mi 

Qr,l 



(K/2 - q 2 y 



Qr,l 



(18) 



Such a form will allow us to easily read off the interference term. 

The matrix element can be expanded in powers of m^/M- for M~ ^> m x , yielding a 
leading (1/M~) 4 behavior. The leading term comes from neglecting the (-K/2 + qi) 2 and 
(K/2 — q 2 ) 2 terms in the propagators compared to M~: 



M 



tree 



l g r/ \ u{qi)T j u PR >L v{q 2 ) 



V2 



Mi 

1R,L 



(19) 



As a check, we can now compute the zero- velocity neutralino annihilation cross section from 
the tree-level XX ~^ QQ9 process, neglecting the final-state quark mass. Squaring the matrix 
element, summing over final-state polarizations and colors and averaging over initial-state 
polarizations, we obtain 



4 



tree 



2gt 



pols 



\9r 



1R 



+ 



\9t\ 

QL 



qi ■ Q2 



(20) 



for a single final-state quark flavor, where we include the contributions from exchange of 
both q L and (Jr. Integrating over the phase space, we find 



Vre\<?(xx -> qqg) 



4 m' 

A 

15 (4?r) 5 



\9r 



QR 



+ 



\gtr 

QL 



(tree level). 



(21) 



This agrees with the result of Drees et al. [5j when their Eq. (2.18) is expanded to extract 
the leading 1/M~ dependence. 

- Q 

The cross section for the similar tree-level process XX ~^ ffl was computed in Ref. [9J] for 



the case of pure photino LSPs and degenerate left- and right-handed squarks, M~ = M- 



QR 



M~. Replacing the color factor with the appropriate electric charge entails a s — > aQ^/A. 
The photino couplings are gg = —g r = y2eQf, so that the couplings |g r /| 4 become Ig^el 4, — ► 

ffj with in j = and 

(22) 



4e Qf. Thus we obtain, for the tree-level (dimension-eight) XX 



Vrcl 



0. 



VreW(XX -> ffl) 



m 



15 M§ 

9 



(tree level). 

Note that this result is a factor of two smaller than that of Ref. . 



C. Interference term 



We now combine the dimension-six and dimension-eight contributions. While the 
dimension-six gluon-splitting amplitude in Eq. (|TT|) contains a purely vectorlike final-state 

9 



quark current, the dimension-eight amplitude contains both vectorlike and axial- vector like 
quark currents: 



9 



x 



Mi 

QR 



+ 



Mi 

1L 



u(qi)T c j v v(q 2 ) + 



Mi 

QR 



Mi 

9l 



u( qi )T c lul5 v(q 2 )}{23) 



Only the vectorlike part of .Mtree will interfere with the dimension-six amplitude. Summing 
over final-state polarizations and colors and averaging over initial-state polarizations, we 
obtain the interference term, 



1 



pols 



MtreeM 



1 split 



\9r 



\9r 



r LL 7 1 r > ' V71 TT" (?1 ' Q-if + (?2 ' q^f 

(24) 

where we have summed over the five light internal quarks q' in -M sp iit of Eq. (fTTj) and over 
the five light external quarks q. Integrating over the phase space, we find the contribution 
to the cross section from the interference term, 



^rcl " int 



a 2 2 2m 2 



327T 3 



X 



\9r 



Mi, 

1r 



m 

1 L 



E 



\9r 



Mi 

IR 



Mi 

1L 



(25) 



We see that the form of this interference term is rather similar to that of the dimension- 
six cross section given in Eq. (|10p. In the special case of degenerate right- and left-handed 
squarks, M~ R 

cross section by a multiplicative constant: 1 



M~ L = M~ we find that the interference term is related to the dimension-six 



^rcl^int — ^relO"dim6 



2 m l x 
3M? 



This result can be compared to the NLO corrections to XX 
flavors and taking the renormalization scale \x = 2m x we obtain 



99 



V Te l(J — WrclOdime 



1 + 



221 (2 



12 



7T 



2 m l x 
3M2 



(26) 

for Nf = 5 light 
(27) 



The second term in the square brackets is the NLO correction from Ref. |8[ and the third 
term is our result for the interference term. Note that for m v ~ 100 GeV the NLO correction 



1 We disagree with the statement of Ref. [2( that the interference is zero for photino-like neutralinos when 
— M—^ . This conclusion was due to a spurious 75 in the loop amplitude matrix element of Ref. Q , 
which resulted in it interfering with the axial-vector part of the dimension-eight amplitude. 



10 



is roughly 60%; for, e.g., M~ ~ 2m x , the interference term cancels off about one quarter of 
the NLO correction. 

Finally, we note that we have not included a term of order a 2 s m^/M~ from the expansion 
of the squark propagator in XX ~* 99 i n powers of 1/M~. This is because the loop integral 
is odd in the argument (m x /M~) with the next term in the series contributing to the cross 
section at order a^m^/M-, which is beyond our present interest. 

III. CONCLUSIONS 

We have computed the cross section for neutralino annihilation XX ~^ QQ9 a ^ order a 2 s /M~ 
arising from interference between the tree-level and loop-induced processes, in the limit of 
zero relative neutralino velocity and ignoring the masses of the five light quark flavors. Our 
result for the total cross section from this interference term is given in Eq. (|25|). We also 
provided the spin-averaged squares of matrix elements for the dimension-six, dimension- 
eight, and interference term contributions to XX ~* QQ9 i n terms of the final-state particle 
momenta in order to facilitate their implementation into Monte Carlos for the computation 
of indirect dark matter detection rates. 
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APPENDIX A 



The matrix element for XX ~^ 99 with on-shell final-state gluons and massless quarks in 
the loop can be expressed as an expansion in powers of m 2 /M~ as follows: 



-Mloop 

where 



a. 



E 



MA 

1r 



ie ^e* a (k 1 )el a (k 2 )k la k 20 , (Al) 



a 2 a 4 



Jin 



F(a) = H 1 1 h 7 TT7 

v ; 9 25 (2n + l) 2 



+ 



(A2) 



with a = m 2 /M~, in agreement with Ref. 
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